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ABSTRACT 

Protons accelerated in the internal shocks of a long duration gamma ray burst 
can escape the fireball as cosmic rays by converting to neutrons. Hadronic inter- 
actions of these neutrons inside a stellar wind bubble created by the progenitor 
star will produce TeV gamma rays via neutral meson decay and synchrotron ra- 
diation by charged pion-decay electrons in the wind magnetic field. Such gamma 
rays should be observable from nearby gamma ray bursts by currently running 
and upcoming ground-based detectors. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mecha- 
nisms: nonthermal 



1. Introduction 



TeV 7-rays have not been yet convincingly detected from a gamma ray burs t (GRB). The 



data from GRB 970417a reported by the Milagrito water Cherenkov det ector (lAtkins et al. 
2000h and from GRB 971110 reported by the GRAND air shower array JPoirier et al-lbooj ) 



lack one of the most crucial pieces of information about those GRBs, namely their redshifts. 
A redshift is necessary to estimate the burst energetics and to properly take into account 
attenuation of TeV 7-rays in the extragalactic background radiation (EBL) fields. Moreover, 
their detection is at the ~ 3a level and hence statistically not very significant. A detection 
at higher significance level has been reported b y the TIBET air showe r array by stacking 
data for a large number of GRB time windows (lAmenomori et al.ll200fl ). Currently, several 
imaging air Cherenkov telescopes (IACTs) such as MAGIC and VERITAS are capable of 
slewing to the GRB direction in the sky prompted by burst alert network. The successor 
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of the recently decomissioned Milagro detector, namely the High Altitude Water Cherenkov 
detector HAWC0, has a high duty cycle and is particularly suitable to detect >TeV 7-rays 
from G RB. Their observa t ions have provided upper limits on 7-ray fluence from several 



GRBs flAtkins et all I2OO5I : [Albert et al]l2006h Mor e powerful detectors covering a wider 



energy range such as AGIS fcrawczvnski et al.ll2007f ) and CTA@ are being planned. 



Theoretical models do not predict TeV 7-rays from the internal shocks of a GRB be- 
cause of a high opacity to e ± pair production with observed keV-MeV energy 7-rays unless 
the L orentz factor of the relativistic bulk motion is very large (IRazzaque. Meszaros fc Zhang 
2004 ) . In the external shock with larger shock radii, TeV 7-rays formed by Compton scatter- 
ing o f shock-accelerated electrons may avoid pair pr oduction (IDermer. Chiang &: Mitman 



20001 ; IWang. Dai &: Lull200ll; IZhang fe Meszaros l200lh . Photohadronic cascades induced by 
shock-accelerated protons ( Bottcher fc Dermer 19981 ) requires intense internal soft photon 



fields that will strongly attenuate TeV 7-rays in the EBL. 

In this Letter we propose a hadronic mechanism to produce TeV 7-rays at the same 
time as the prompt keV-M eV emission. If protons are accelerated in the internal shocks 
of a GRB (jWaxmanl I1995T ). they are expected to interact with observed keV-MeV pho- 
tons to produce neutrons (th — » mr + ) which may escape the shock region a s cosmic- 
rays (jWaxman fc Bahcall 1997 ; Rachen fc Meszaros 1998 ; Dermer fc Atoyan 2003 ). A frac- 



tion of these cosmic rays will interact with particles in a surrounding dense ste llar wind as the 
proge nitor star is expected to undergo substantial mass loss before explosion (I Chevalier &; Li 
19991 ). Neutral pions from secondary nuclear production promptly decay to produce very 
high energy 7-rays, while charged pions decay to produce electrons which emit synchrotron 
radiation in the magnetic field of the stellar wind. These 7-rays at TeV energy could be 
detected if they avoid substantial absorption in the source environment as well as 

We find that TeV 7-rays formed by nuclear interactions of escaping neutrons with stel- 
lar wind particles may be detected from a nearby GRB by current and upcoming 7-ray 
Cherenkov telescopes. 



2. GRB internal shocks and cosmic-ray escape 

The GRB internal shocks take place over a wide range of fireball radii, depending 
on the 7-ray variability time scale t v ~ 10~ 3 s and the Lorentz factor of the bulk out- 
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flow Tfe > 10 2 5 r( )) 2.5 for typical long GRBs. For observed non-thermal emission, the radii 
where the internal shocks occur need to be larger than the jet photospheric radius r p h = 
(crT-^ 7 ,iso)/(4vre e r^m p c 3 ) ~ 3.7 x 10 11 L 7j51 e~i i r^2 5 cm, at which the fireball becomes opti- 
cally thin to Thomson scattering. Here we use an isotropic-equivalent 7-ray luminosity of 
-^7,iso = 10 51 L 7i 5i ergs s -1 and a kinetic luminosity of Lk.iso = Lj^e' 1 , where e e = 0.1e ej _i 
is the fraction of kinetic energy converted to 7-rays (assuming a fast-cooling scenario). 
With these parameters we calculate a pre-shock electron and baryon number density of 
n e — n 'p — -^7,iso/(47re e rg h r 2 m p c 3 ) in the comoving frame. We denote the variables in the 
comoving plasma (observer's) frame with (without) primes. 

For our modeling purpose, we assume a shock radius of r s h = 10 12 r s h 12 cm. The turbu- 
lent magnetic field strength in the shock region, assuming that the magnetic energy density 
B'JJ8ti is a fraction e# = 0.1e_e,-i of the fireball's kinetic energy density n' p m p c 2 , is 

B' h » 8.2 x 10 5 ^L^e^Sr^T^ G. (1) 

The protons and electrons are assumed to be accelerated via a Fermi mechanism by this 
magnetic field. 

Characteristic synchrotron photons radiated by the population of electrons with a min- 
imum Lorentz factor j' emin — e e {m p /m e ){T re \ — 1) is one of the leading models to produce 
observed 7-rays. With the parameters adopted here and a relative Lorentz factor between 
two colliding shells r rc i rj 3, the observed characteristic synchrotron photon energy is 

,'2 



600(1 + z)-^J*e%l x Llj*r£ l2 keV. (2) 



(3/2)(B'jB Q ) le 2 min T b m e c 2 /(l + z) 

-1/2 -1 

v 7,51 r sh,12 

Here Bq = m 2 c^ /qh ~ 4.414 x 10 13 G. The observed photon energy at the peak of the 
eF £ spectrum, e p k, varies from burst to burst, however, there exist several phen omenological 



relations connecting the peak photon energy to other burst para meters (e.g., lAmati et al. 



20021 ; iGhirlanda. Ghisellini fc Lazzatill2004l ; IWillingale et al.ll2008l ). Here we adopt a relation 



between the peak 7-ray luminosity and a ch aracteristic photon energy, which in turn is related 



to e p k as found by IWillingale et al.l (120081 ) . We rewrite this relationship as 



e pk ^ 650(1 + z)- 1 L° 1 % keV. (3) 
Note that this is close to the value of the synchrotron photon energy in equation (j2J). 

Following the phenomenological broken power-law fits, we write the comoving photon 
spectrum as 



n' 7 (e') 9* n'^ pk T b /[e pk (l + z)] 
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f «a/4 k )~ a (^/4a) 3/2 ; ^ < 4 

x c^/ey-"; <a < ^ 4k ( 4 ) 

I ( £ / £ pk) ) £ max > £ ' > £ p k 5 

where (e^, e^ax) = (10~ 2 ' 5 , 10 6 ) keV are respectively the synchrotron self-absorption and 
maximum photon energies. The fitted values for the power-law indices are (a, (3) = (1, 2.3). 
We calculate the peak photon number density L ltiso / (^nr^cTbe^) , including a bolometric 
factor of ~ 3 and using equation as 

< pk « 2.7 x 10 18 (1 + z)L°^T^ 5 r; h % cm" 3 . (5) 



The energy gained by the protons is proportional to the time t pacc ~ <f)E' p / qB' sh c, they 
spend in the shock region. The maximum energy is typically obtained by requiring that 
this time with <fi > 1 to be equivalent to the smaller of the fireball expansion or dynamic 
time t^ yn ~ r sil /2cT b and the energy loss time scale t ploss . The synchroton energy loss time 
scale t' p syn = 6irm p c 3 / {aTm 2 E' p B'^) is typically the shortest for internal shocks. A maximum 
cosmic ray proton energy can thus be obtained as 

e ~7xin 10 f 1/4 r 3/2 r 1/2 r 1/4 r 1/4 fi P v m 

-frp.max ~ ' X iU e e,-l L b,2.5 T sh,12 6 B -1^7,51 ^ e V , (0 ) 

for t^ cc = t' syn . The differential spectrum (e.g., in units of cm -2 s" 1 GeV -1 ) of cosmic ray 
protons, if they could escape freely from the fireball at a luminosity distance di, may be 
written as 

J P (E P ) = L 7iiso /[47r4e e £ p 2 \n(E m /T b m p c% (7) 

where we have assumed a typical N(E) oc E~ 2 spectrum generated in a mildly relativistic 
shock. 

Shock-accelerated protons are expected to be confined in the GRB fireball by the mag- 
netic field. Particles can, however, escape directly from the internal shock region when 
protons convert to neutrons through — ■> n7c + interaction. The rate of scattering with 
observed 7-rays, assumed to be isotropically distributed in the GRB fireball, is given by 

KM) = A / defaM) / , ds'^l. (8) 

Here e' r = 7 p e'(l — /3 p cos9) is the photon energy evaluated in the proton's rest frame for 
the angle 6 between the directions of the energetic proton and target photon, and e' th = 
m n c 2 + m 2 c 2 /2m p is the threshold photon energy for pion production. The dominant neutron 
production channel is py — > nii + with an intermediate A(1232) baryonic resonance produc- 



tion. The cross section formula may be written in the Breit-Wigner form (IMiicke et al.ll2000l ) 
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2 \21— 1 



and the peak cross section is given by <7 p k = 4.12 x 10 
nance width is Ta = 0.11 GeV. 



where s = mic A +2e' r m p c 2 , <To 



-28 cm 2 £ i 



r,pk 



3.11 xlO" 29 cm 2 



0.3 GeV. The reso- 



For the photon spectrum in equation (j3J), and taking an approximate cross section given 
by a pl {e' r ) = cr p k for e' k < e' r < e' k + Ta, the scattering rate in equation ([8]) simplifies to 



ccr p kn 7ipk 



(/3 2 -l)£ 



'/3-1 
r.pk 
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(9) 



2 l/2 £ '5/2 



2-v' 2 r' 2 



57; 



'5/2 '5/2 ) /p 



V > £a 

"t> ^ f' 



for a = 1. Here e' A = (e' r pk — T/^rripC 2 . The rate K{^' p ) times the dynamical time t' dyn ~ 
r s h/2cFb represents the opacity 7^(7') for jry scattering. For target photons with energies 
below e' pk , the p7 scattering rate of high-energy protons is almost constant with energy when 
a — 1. Therefore T P7 (7p) ~ -^(7 p )^d y n ~ 1 in this plateau region of 7' where ^(7^) refers to 
the accurate rate, equation (|S]), or the approximate rate, equation respectively. 

The flux J n (E n ) of neutrons from jry interactions depends crucially on the shocked- 
fireball radius, as discussed above. Assuming that they don't interact further in the fireball, 



J n (E n ) = J p (E n /y)t' dyn K(E n /yT b )/2y 



(10) 



where y = E n /E p « 0.8 is the fraction of proton energy given to a secondary neutron, where 
a mean inelasticity is used. The injected primary proton flux (topmost thick dashed line) and 
escaping neutron flux are plotted for the full numerical calculation (thin solid line overlayed 
on dots) and for the approximate expression (dots) in Fig. [TJ 



3. Cosmic ray interaction with stellar wind 



Massive stars such as GRB progenitors lose mass by blowing a spherically symmetric 
and steady wind. We assume a nominal mass loss rate of M w = 10~ 4 ' 5 M U) _4. 5 M yr -1 and 
a wind velocity of v w = 10 8 f Wj8 cm s_1 - The volume density of particles in the wind is 
M w / (4irr 2 v w m p ) and the column density at a radius r = r sh = 10 12 r 12 cm is 



9.5 x 10 23 M 



-4.5 Vs^ 1 cm 2 - 



The stellar wind may have high magnetic field, as has been suggested by many authors (e.g., 
Volk fc Biermannl Il988l ; iBiermann fc CassineUil Il993l ) . We assume for simplicity that this 
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field is in equipartition with the wind kinetic luminosity MyjV^/2 (jWang et al.ll20071 ). so that 

B w » 141 < 2 _ 1 M^! 4 . 5 ^rr 2 1 G , (12) 
where wb = O.Iwb -i is the equipartition parameter. 



Cosmic ray neutrons escaping from the GRB internal shocks can interact with dense stel- 
lar wind particles and produce secondary pions, kaons, and higher-order resonances through 
pn interactions. The neutron decay radius is 3> cr^T b 10 16 I\ 2 .5 cm, so that neutrons that 
do interact with wind particles do so before they decay, and will therefore make beamed 
secondaries that would be directed along the GRB jet. Neutral pion and eta mesons decay 
almost instantaneously to produce ultrahigh energy 7-rays. 

The 7-ray flux from pn interactions of neutrons with stellar wind can be calculated from 
the expression 



JJE~) 



f 1 dx 

1 <Jn 






E~ 


Jo x 


X 


X 



(13) 



Here a P p(E p ) is the inelastic pp cross-section, x = E^/E n is the fractional 7-ray energy and 
F 7 (x;i? 7 ) is the 7-ray yie ld function from neutral meson decay s. We use the Y 7 (x;E 7 ) as 



recently parametrized by iKelner. Aharonian Bugayovl (120061 ) of the SIBYLL code which 
include 7-ray production from both tt° and rf decays. Note that the charged lepton flux from 
pion decays may also be calculated using equation (fT3j) with a change of subscript 7 — ¥ e and 
using the appropriate yield function. The 7-ray (thin dash-dotted line) and electron (thin 
dashed line) source fluxes are plotted in Fig. CE] 

Comparing the synchrotron cooling time scale t e ,s yn = (3/2)/?, 2 (_B^/_BQ) _2 (r e m e c£ , e )~ 1 
for 7r ± decay in the wind magnetic field given by equation (II 2p . with the observed 7-ray 
variability time t v ~ r s h/2r 2 c, we find that electrons with -E e ,min > 8 x 10 4 GeV radiate 
away a large fraction of their energy. Here r e is the classical electron radius, and we let the 
perpendicular magnetic field B± = B w . The total synchrotron power emitted by an electron 
is given by P = (2/3)(r e /h 2 )(B^/BQ) 2 E^m e c with a characteristic photon energy E c = 
(3/2)(B±/BQ)E 2 /m e c 2 , similar to the expression in equation (j2J). To a good approximation 
we can assume that the total power ATtd 2 L t v PE e J e (E e ) is emitted in photons of energy E c . 
The corresponding synchrotron flux by the electrons is therefore given by 



E^Jsyn^E^) 



B„ 



(3/2m e ) 5 /2£2 y Bq ) J e(£)- 



1/2 



(14) 



Here £ = [(2/3)(B w / ' Bo)E 1 m e c 2 ] 1 / 2 . Note that the synchrotron energy loss formula assumed 
here applies in the classical limit defined by the parameter \ — (3/2) Bo)(E e /m e c 2 ) <C 1. 
Thus the minimum and maximum synchrotron photon energies for the parameters adopted 
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here may be calculated as E 1 



{21/2){B Q /B w fh^ b /[mlc 



2Jl„2 



y 7,min 

_ 2 GeV, assuming \ 
plotted in Fig. [T] with an exponential cutoff above -E 7imax 



shJ 



60 GeV and -E 7 , m ax — 



(2/3) (Bg/BJxW 



10 X-2- The synchrotron flux is 



Compton losses on the scattered stellar radiation field can be shown to be small com- 
pared with synchrotron losses. The energy density of scattered stellar photons from the 
pre-burst star is ~ L ls r w /A'Kr 2 c = 260 L^r w /r\ 2 ergs cm~ 3 , where = 10 38 L* j3 g ergs s -1 
is the pre-burst stellar luminosity and r w is the Thomson depth of the wind. This is smaller 
than the magnetic field energy density B^/871 given from the expression for B w in equation 
(|12p . even for a luminous pre-burst star. Klein-Nishina effects will make the Compton losses 
even smaller. TeV 7 rays might also produce pairs with the stellar photons through 77 
interactions. The optical depth of TeV photons to 77 pair production with stellar photons 
with mean energy e* can be written as r 77 ~ ((7x/3)n p ?i(e*)r ~ 4 x 10 _5 L* i 38TT/[ri2(e*/eV)]. 
As can be seen, this process can be neglected. 

High energy 7-rays are also subject to absorption with photons of the EBL while prop- 
agating fro m the source to Earth. The opacity r 77 ~ 1 for « 600 GeV photons from a source 
at z « 0.1 (iRazzaque. Dermer fc Finkd 120081 ) . To calculate the opacity, we assumed that 
the background radiation field consists of three components — cosmic microwave background, 
infrared a nd optical photo ns — represe nted by a blackbody spectrum, a modified blackbody 
spectrum (jDermerl l2007bl ) and a fit (IRazzaque. Dermer fc Finkd 120081 ) . respectively. The 
final emerging 7-ray spectrum (thick solid curve) is plotted in Fig. [TJ The absorbed 7-rays 
can induce a pair cascade and give rise to a long duration com ponent after the burst is over if 
the in tergalactic magnetic field is sufficiently weak, < 10~ 16 G ( IRazzaque. Meszaros Zhang 
20041 ) . 



Results and Discussion 



Figure [JJ shows the results of our study. We assumed a GRB luminosity distance of 
d,L = 455 Mpc [z = 0.1), with all scaling parameters equal to unity (most importantly, r 12 = 1 
and M w -4.5 = 1). The "observed" 7-ray spectrum is calculated from the "e synchrotron" 
and "7T°, 77 — > 7" components after taking into accou nt absorption in the EBL. Also shown 
in Fig. [1] are the detection sensitivites of the MAGIC ([Albert et al.ll2006l ; IScapin et al.ll2006l ) 
and HAWC (see footnote 4) detectors. For MAGIC, we used their 60 s 5cr GRB sensitivity 
of 5.8 Crab between 80 GeV - 350 GeV and 1.8 Crab between 350 GeV - 1 TeV. For HAWC 
we used their 10 s 5cr GRB sensitivity within 0-10 degrees of the azimuth. 



As shown in Fig. (TJ a typical long duration GRB inside a stellar wind environment may 



- 8- 



be detected by IACTs with rapid slewing capability such as MAGIC, VERITAS or HESS 
and by the upcoming HAWC detector. If all long duration GRBs within z ~ 0.1 have dense 
stellar wind as modeled here, then the expected detection rate in upcoming TeV detector 
would be > 1 burst yr" 1 , using a GRB rate of 2 Gpc~ 3 yr _1 . HAWC would be sensitive for 
GRBs with z < 0.2, when distance and EBL effects become important. 

In the context of the internal shock model, a prompt TeV emission signal will result 
from the first pair of colliding shells that form neutrons which escape and then interact with 
particles in the wind. As the merged shell moves out along the GRB jet, it sweeps up wind 
material, so that subsequent escaping neutrons will no longer have target wind particles with 
which to interact. The column density of material does not however change, so that neutrons 
formed by further pairs of colliding shells still have a significant target column density with 
which to interact and make TeV radiation. Thus the duration of the prompt TeV signal in 
this model corresponds to the duration of the prompt phase associated with colliding shells. 

Predicted TeV 7-ray emiss ion in the early afterglow phase, either by ha dronic interac- 

tions ( Bottcher fc Dermerlll998) or synchrotron self-Co mpton (SSC) emission ( Dermer. Chiang fc Mitman 



2000 ; IWang. Dai fc Lull200ll ; IZhang fc Meszarod 120011 ). is expected to last much longer than 



the prompt TeV emission considered here. A leptonic SSC origin of TeV radiation formed 
by an external shock will correlate with the lower energy synchrotron radiation with a peak 
energy that becomes smaller as the blast wave decelerates. By contrast, the TeV emission 
formed by the process considered here will end when the central engine becomes inactive. 
The TeV 7-ray flux predicted in this work should correlate with the activity of the central 
engine as reflected by the MeV emission from a GRB, though delays could arise from the time 
required to accelerate protons to ultra-high energies. High-energies neutrinos are formed di- 
rectly from photopion-producing interactions in the internal shocks (in the TeV-PeV energy 
range), as well as from pn interactions (in the PeV-EeV energy range) in the stellar wind. 
Joint detection of prompt high-energy neutrinos and prompt TeV radiation would provide a 
new method to probe the environment in the vicinity of GRBs. 



The work of S.R and C.D.D. is supported by the Office of Naval Research. 
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Fig. 1. — Spectra of shock-accelerated cosmic rays, neutrons escaping from the GRB fireball 
(thin solid line shows the numerical calculation using the full jry — > A cross section formula, 
and the dots show equation Q) using a constant cross section), n° and rf decay 7-rays 
and ii^ decay electrons. The e synchrotron radiation is calculated from pion-decay electrons 
in the wind magnetic field, and the thick solid line is the emerging 7-ray spectrum after 
absorption in background radiation fields. Also plotted are the MAGIC and HAWC detector 
sensitivities. 



